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CN ■ Abstract 



We are interested in the evolving genealogy of a birth and death process with trait 

structure and ecological interactions. Traits are hereditarily transmitted from a parent to its 

offspring unless a mutation occurs. The dynamics may depend on the trait of the ancestors 

rvq I and on its past and allows interactions between individuals through their lineages. We define 

CN ■ an interacting historical particle process describing the genealogies of the living individuals; 

it takes values in the space of point measures on an infinite dimensional cadlag path space. 

This individual-based process can be approximated by a nonlinear historical superprocess, 

under the assumptions of large populations, small individuals and allometric demographies. 

Because of the interactions, the branching property fails and we use martingale problems and 

■■M I fine couplings between our population and independent branching particles. Our convergence 

Cd ■ theorem is illustrated by two examples of current interest in biology. The first one relates 

the biodiversity history of a population and its phylogeny, while the second treats a spatial 

model with competition between individuals through their past trajectories. 

^J . Keywords: Nonlinear historical superprocess; Genealogical interacting particle system; Limit 

^S| i theorem; Evolution models. 

^n ■ AMS Subject Classification: 60J80, 60J68, 60K35. 

(N 

!>. 

5 ■ 1 Introduction 

The evolution of genealogies in population dynamics is a major problem, which motivated an 
abundant literature and has applications to evolution and population genetics. Our purpose 
K> ! here is to generalize the existing models by emphasizing the ecological interactions, namely the 

C^ ' competition between individuals for limited resources. In this paper, we construct a structured 

birth and death process with mutation and competition, whose dynamics depends on the past. 
Each individual is characterized by a vector trait x € M*^, which remains constant during the 
individual's life and is transmitted hereditarily unless a mutation occurs. The birth and death 
rates of an individual can depend on the traits of its ancestors or on the trait's age in the lineage. 
Moreover individuals interact with each other. 

Here, we are interested in keeping track of the genealogies of individuals with small weights, 
in large populations with allometric demographies (i.e. short individual lives and reproduction 
times). At a given time t, we associate to each individual its lineage until t defined as a function 
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which associates to each s <t the trait of its ancestor living at s (see Fig. [T|) . Since each individ- 
ual keeps a constant trait during its hfe, these hneages are cadlag piecewise constant paths. The 
population is described by a point measure on the lineage space. To reflect allometric demogra- 
phies, we introduce a parameter n which scales the size of the population. The individuals are 
weighted by 1/n to keep the total biomass of constant order when n varies. Moreover, lifetimes 
and gestation lengths are proportional to this weight. Hence birth and death rates are of order n 
and preserve the demographic balance. Mutation steps are assumed to have a variance of order 
1/n. The trait evolution is driven by mutations and competition between individuals. 

traits 

• 

• c 
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time 

Figure 1: Example of a path constituted with ancestral traits. 

We study the convergence of these processes in large populations, when n tends to infinity. 
We proceed with tightness-uniqueness arguments inspired by Dawson and Perkins [TS] for the 
historical super-Brownian process (without interaction). In our model, there is no branching 
property, impeding the use of Laplace's transform arguments. We introduce a new infinite di- 
mensional martingale problem and use fine comparisons between our population and independent 
branching particles. Moreover, the measure-valued processes on path space that we study are 
discontinuous, which makes it necessary to define a new class of test functions. In the limit, we 
obtain a nonlinear historical superprocess, generalizing the work of Perkins |36j . In particular, 
past dependence and competition can be taken into account in the lineages. Examples are given, 
suggesting that the historical processes may open the way to models of evolution without the 
assumption of rare mutations and time scale separations. Let us remark that our model allows 
the description of both genealogies and population densities in the forward physical time, with 
variable population size including extinction phenomena. This can be a first step to reconstruct 
the past biodiversity from the phylogenies of living species with ecological interactions. 

From this convergence result, we study the dynamics of the distribution of genealogies. We 
show that it can be seen as the dynamics of a diffusive particle system resampled in such a way 
that individuals with larger allometric (diffusive) coefficient or with higher growth rate are more 
likely to be chosen. A particular case is when these two rates are constant. Then, the intensity 
of the distribution of the genealogies corresponds to a Fleming- Viot process with genealogies, 
similar to the one introduced in Greven-Pfaffelhuber- Winter |24j . In the case without interaction 
or with logistic competition without past or trait dependence, one can also obtain Feynman- 
Kac's formulas for the mass process renormalized by its expectation. This formula can help 
in getting easier simulation schemes for certain demographic quantities. Let us note that more 
generally, in absence of interaction, coalescent processes provide another modeling of genealogies 
(see Berestycki [5] and references therein for a survey). Depperschmidt, Greven, Pfaffelhuber 
and Winter [16^ [M] represent genealogical trees (for instance A-coalescents) as marked ultra- 
metric spaces and the absence of interaction allow them to study genealogical distances by using 
Laplace transforms. In a recent work. Barton, Etheridge and Veber [3j study the genealogies of 
a spatial version of A-Fleming-Viot processes and of their various limits. All these models allow 
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the incorporation of selection and mutation (see also e.g. |3H [3l [20]) but not of competition 
between individuals, which is the main aim of our work. 

In Section 2, we construct the historical particle system whose dynamics depends on the past. 
The diffusive limit in large population is obtained in Section 3 under large population and 
allometric demography assumptions. We prove tightness of the laws of the historical particle 
processes with new arguments and test functions accounting for interactions and jumps. We then 
identify the limiting values as solutions of a nonlinear martingale problem for which uniqueness 
is stated. The distribution of genealogies is studied in Section 4 and examples are carried in 
Section 5. The first example deals with an evolution model of adaptive dynamics with local (see 
|17j ) competition. The historical superprocess dynamics shows that for a range of parameters, 
the population separates into groups concentrated around some trait values. This cannot be 
observed by the unique information on the (classical) superprocess as noticed in Figure [2j This 
example suggests that the historical processes may allow one to understand evolution without the 
assumption of rare mutations and time scale separation. For example, our model could be a basis 
to reconstruct the past biodiversity from the actual molecular phylogenies. Until now the models 
which have been used do not take into account the interaction between individuals (see Morion 
et al. [S])- The second example is a spatial model (see [1] or [361 P-50]): particles consume 
resources where they are locating. Therefore, the offspring arriving in previously habited regions 
are penalized. We will see that this tends to separate the cloud of particles in several distinct 
families whose common ancestor is very old. 

Notation: For a given metric space E, we denote by B^; = D(M_(_,£') the space of cadlag 
functions from M4. to E. For ii^ = M, we will use the more simple notation D = D(R_|_,M). 
These spaces are embedded with the Skorohod topology which makes them Polish spaces, {e.g. 
[elEHlES], see also (|Bl]) in appendix). 

For a function x € B^; and t > 0, we denote by x* the stopped function defined by x*(s) = x(s At) 
and by x*~ the function defined by x*~(s) = liiUr^t x''(s). We will also often write xt = x{t) for 
the value of the function at time t. For y,w G Oe and t £ M+, we denote by (yltlw) € B^; the 
following path: 

(^1^1") = 1 u;„_. if n>t. (^-1) 

When the path w is constant with Vu € M4., Wu = x, we will write (y|t|x) with a notational 
abuse. 

We denote by M.f{E) (resp. M.'^{E), 'PiE)) the set of finite measures on E (resp. of point 
measures renormalized by 1/n, of probability measures). These spaces are embedded with the 
topology of weak convergence. 

2 The historical particle system 

In this section, we construct the finite interacting historical particle system that we will study. 
Trait-structured particle systems without dependence on the past have been considered in 
Fournier and Meleard [23] or Champagnat et al. [10]. For populations with age-structure, 
we refer to Jagers [271 ES] and Meleard and Tran |33l [39] for instance. Here, we are inspired by 
these works and propose a birth and death particle system where the lineage of each particle, 
i.e. the traits of its ancestors, is encoded into a path of BjRd. 



2.1 Lineage 

We consider a discrete population in continuous time where the individuals reproduce asexually 
and die with rates that depend on a hereditary trait and on their past. Each individual is asso- 
ciated with a quantitative trait transmitted from its parent except when a mutation occurs. The 
rates may express through the traits carried by the ancestors of the individual. One purpose 
is for example to take into account the accumulation of beneficial and deleterious mutations 
through generations. 

Individuals are characterized by a trait a; G M'^. The lineage or past history of an individual is 
defined by the succession of ancestral traits with their appearance times and by the succession 
of ancestral reproduction times (birth of new individuals) . To an individual of trait x born at 
time Sm, having m — 1 ancestors born after t = at times 5"! = < ^2 < • • • < Sm-i, with 
Sm~i < Sm, and of traits (xi, X2, • • • Xm-i), we associate the path 

m—l 

yt=Yl ^As,<t<s,+-, + xls^Kt- (2.1) 

This path is called the lineage of the individual. We denote by C the set of possible lineages of 
the form ()2.ip . Since a path in C is entirely characterized by the integer m and the sequence 
(0, xi, . . . , Sm-i,Xm-i, Sm, x) of jump times and traits, it is possible to describe each element of 
£ by an element of N x |Jj^gj^(M+ x M*^)™ embedded with a natural lexicographical order. 

2.2 Population dynamics 

Let us introduce a parameter n G N* = {1,2,...} scaling the carrying capacity, when the total 
amount of resources is fixed. To keep the total biomass constant, individuals are attributed a 
weight 1/n. The population is represented by a point measure as follows: 

^t ■■= IJ2^\^ ^ ^P(^^ ^ M},{n^,), (2.2) 

where NJ^ = n {X^, 1) is the number of individuals alive at time t. 

The reproduction is asexual and the offsprings inherit the trait of the ancestor except when a 
mutation occurs. Death can be due to the background of each individual or to competition 
with the other individuals. We consider allometric demographies where lifetimes and gestation 
lengths of individuals are proportional to the biomass. Thus, birth and death rates are of order 
n, but respect the constraint of preservation of the demographic balance. Also, the mutation 
steps are rescaled by 1/n. 

Let us now define the population dynamics. For n € N*, we consider an individual characterized 
at time t by the lineage y G Djgd in a population X^ G D;v4"(d d)- 

Reproduction: The birth rate at time t is b"'{t,y), where 

6"(t,y) = nr(t,y) + 6(t,y). 

The function 6 is a continuous nonnegative real function on M-|_ x Bjgd. For instance, b can be 
chosen in the form: 

b{t,y)=B(^J yt-sMds)^ (2.3) 
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where i? is a continuous function bounded by B and ff, is a Radon measure on M+. The 
allometric function r is given by: 



rit,y) = R 



( f yt-sMds)j (2.4) 



where R is continuous and bounded below and above by ^ > and R > 0, and where Ur is a 
Radon measure on M-(_. We also assume that R is Lipschitz continuous, implying that i?^'^ is 
also Lipschitz. 

When an individual with trait yt_ gives birth at time t, the new offspring is either a mutant 
or a clone. With probability 1 — p G [0,1], the new individual is a clone of its parent, with 
same trait yt_ and same lineage y. With probability p G [0,1], the offspring is a mutant of 
trait yt_ + h, where h is drawn in the distribution k'^{h) dh. To this mutant is associated the 
lineage {y\t\yt_ + h). For the sake of simplicity, the mutation density k'^{h) is assumed to be 
a Gaussian density with mean and covariance matrix o"^ Id/n. However the model could be 
generalized for instance to mutation densities k"'{yt_, h) with dependence on the parent's trait. 
Let us introduce the notation: 

K^'idh) = pk''{h)dh + (1 - p)6o{dh). (2.5) 

Example 2.1. (i) If tabids) = 6o{ds), then Jq yt-s^b{ds) = yt is the trait of the individual of the 
lineage y living at time t. 

(ii) If I'bids) = e~"'^ds, with a > 0, then f^ yt-s'-'bids) = /q e^^'^^^'^'ygds. This means that the 
traits of recent ancestors have a higher contribution in the birth rate of the individual alive at 
time t. Such rates may be useful to model social interactions, for instance cooperative breeding 
where the ancestors contribute to protect and raise their descendants. When ancestors have 
advantageous traits, they may help their offspring to reproduce in more favorable conditions 
and increase their birth rates. □ 

Death: To define the death rate, let us consider a bounded continuous interaction kernel 
U G C;,(M+ X D^^,]R), a bounded continuous function D on M+ x H^d and a Radon measure z/^ 
weighting the influence of the past population on the present individual y at time t. The death 
rate is 

d^{t,y,X^) = nr{t,y)+d{t,y,X^), 

where for a process X G Omf(b> a)^ 

d{t,y,X)=D{t,y)+ f f U{t,y,y')Xt^s{dy')Mds). (2.6) 

The first term with function r allows us to preserve the demographic balance. The term D{t, y) 
is the natural death rate, while U{t, y, y') represents the competition exerted at time t on the 
individual of lineage y, by an individual of lineage y' alive in its past. We assume that: 

3D > 0, Vy G B, Vt G M+, < D{t, y) < D, 

3U, U>0, Vy, y' G B, Vt G E+, < U{t, y, y') < U. (2.7) 

Example 2.2. Model with asymmetrical competition: 

If we choose Z^ = 0, Vd{ds) = 5o{ds) and the asymmetric competition kernel proposed by Kisdi 

m 

2 / 1 
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''^'^y^y'~^ = K{'-TT—^^^ ^2-^ 



then, the death rate becomes: 

Remark 2.3. A generahzation to physical age structure as considered in [32] is possible, pro- 
vided we extend the trait space by giving a color to each individual. The colors are independent 
uniform [0, l]-valued random variables drawn at each birth. The lineage of colors of an individ- 
ual is a path c S D(M+, [0, 1]) and allows the definition of the birth date of an individual alive 
at time t, 

Tc,t = inf{s < t, Cs = ct} = sup{s < t, Cs ^ cj (2-10) 

The age of the latter individual at time t is given by 

a{t) := t - Tc,t. (2.11) 

It is a cadlag function, discontinuous at the birth times. □ 

2.3 Construction of the historical particle process 

Following the work of Fournier and Meleard |23j . the population process is obtained as a solution 
of a stochastic differential equation (SDE) driven by a Poisson Point Measure (PPM), describing 
the dynamics of (Xf )igiR^, for any n € N*. The measure representing the population evolves 
through the occurrences of births and deaths. Since the rates may vary with time, acceptance- 
rejection techniques are used to obtain these events' occurrences by mean of PPMs. According 
to birth and death events, Dirac masses are added or removed. This construction provides an 
exact simulation algorithm that is extensively used in Section | 




Remark 2.4. As in Fournier and Meleard [23], the map Y = (y*)jgN. from IJnGN* -^p(^) in 
C is defined for n and N in N* by 

otherwise, 

where the individuals are sorted by the lexicographical order. This will be useful to extract a 
particular individual from the population. When there is no ambiguity, we will write Y^ instead 
of Y'^{X) for a point measure X. 



Definition 2.5. Let us consider on the probability space (il, J-", P): 

1. a random variable Xq € A4^{B)^d) satisfying E {{Xq, 1)) < +oo and such that the support 
of Xq contains a.s. only constant functions, 

2. a PPM Q{ds, di, dh, dO) on M+ x £; := M+ x N* x M^ x R+ with intensity measure ds (g) 
n{di) ® dh® dO and independent from Xq, n{di) being the counting measure on N* and 
ds, dh and dO the Lebesgue measures on R+, M'' and M+ respectively. 

We denote by {J-^)teM.+ the canonical filtration associated with Xq and Q, and consider the 
following SDE with values in M^ip^d): 



X? = XR+- t I 1 



'* ~ ^ n j /_-^{i<"(-^?_.i)} 



\Y^\s\Y^_+h)%<m^(i,s,h) + ^Y^'^m'l{i,s,h)<e<m':^{i,s,h) 



^Yi'^mr^{i,s,h)<e<mr^{i,s,h,X''' 



Q{ds,di,dh,de), (2.12) 
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where 



m'lii,s,h) =pU'{s,Y\X^_))k'^{h), 

m^{i, s, h) = m1{i, s, h) + {l- p) 6"(s, Y\X^_ ))r (/i) 

m^{i, s, h, X"'^- ) = m^{i, s, h) + (f{s, Y\X^_),X''''- )k''{h). 



(2.13) 
D 



Existence and uniqueness of the solution (X")tgR^ of SDE ()2.12p are obtained for every n G N*, 
from a direct adaptation of [23], as well as the next proposition concerning moment estimates 
and martingale properties. 



Proposition 2.6. Let us assume that 



supE((Xo",l)3) < + 

neN* 



OO. 



(2.14) 



Then, 

(i) For all T > 0, 



sup E( sup {XJ",!)^) < +00. 



(2.15) 



(ii) For a bounded and measurable function if, 



(Xr, v) = {XS, ^)+Mr + ds X^{dy) 



nr{s,y)[ / ip{y\s\ys + h)K'"-{ys,dh) - ip{y) 

+b{s,y) f ip{y\s\ys + h)K^iy,,dh)-d{s,y,iX^rMy) 

JR'I- 

where M"^'^ is a square integrable martingale starting from with quadratic variation 



nJo 



(2.16) 



{nr{s,y) + b{s,y)) / ip'^{y\s\ys + h)K''{ys,dh) 

+ {nr{s,y) + d{s,y,{X^r))^\y)]x^{dy)ds. (2.17) 



2.4 Examples 

In this section, we give two examples of applications of historical processes in biology. 

2.4.1 A model ^vith competition for resources 

We first introduce a model of adaptation with competition for resources that has been considered 
by Roughgarden [37j, Dieckmann and Doebeli [17J, Champagnat and Meleard [12j. In this model, 
the trait x € [0, 4] can be thought as being a (beak) size. The birth and death rates are chosen 
as 

r(t,y) = l, 6(t,y) = expf-%^), d(t,y,X)= [ exp ( - ^y'~ f' )xidy'). 



2al ;' 



2a 



u 



Here, there is no dependence on the past and 7(5, y,X*) = 7(7/3, X^). Non-historical superpro- 
cess renormalizations in such models have been investigated for instance by Champagnat et al. 
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|llj . The birth rate is maximal at x* = 2 and there is a local competition with neighbors of 
close traits. A flat competition kernel {au = +00) will make evolution favoring individuals with 
maximal growth rate x*. For au < +00, Champagnat and Meleard [12] proved that, depending 
whether ab < crjj or at > crjj, either evolution concentrates on individuals with trait 2 or favours 
several groups concentrated around different trait values. 

(a) 




(b) 





Figure 2: Dieckmann-Doebeli model, a = 0.4, ab = 0.4, xq = A, p = 0.5, n = 300. The 300 particles 
are started with the trait 1.5. (a) au = 0.3 ; (b) ajj = 0.45. When a^ > ajj , we observe in the historical 
superprocess a separation of the population into subgroups concentrated around different trait values and 
attached to the MRCA by a single branch. This genealogical separation 0/ the population into a small 
number of families (2 in (a)) cannot be seen on the superprocesses (left column). When a^ < ajj , then 
the families are less distinct and the population has a shorter persistence. 



In Figure [21 we represent the dynamics of the population process including the extinct branches 
and the historical particle process. These two representations complement each other. The 
population process gives the evolution of the number density of individuals, while the historical 
particle process reveals the evolving phylogenies and provides the ancestral paths of living par- 
ticles. The simulations are obtained for a large n and show two different behaviors. It is known 
that the limit of the population process is a superprocess (left pictures). Our aim is to under- 
stand the limit of the historical particle process (right pictures), i.e. the nonlinear historical 
superprocess. 

In the superprocess, the individual dimension is lost and ancestral paths can not be read from 
the sole information of the support of the measure. The latter represents the dynamics of the 
population's biodiversity. The interest of the historical superprocess is to provide genealogical 



information. Indeed, one can stress that the historical process at time t restricted to paths up 
to time t — e has finite support, as observed in Figure [21 This result is known for the historical 
super Brownian motion |151 Section 3]. 

In the simulations of Figure [21 we see that when ab > (Tjj, two families appear, in the 
sense that in the genealogical tree (on the right) the most recent common ancestor (MRCA) is 
separated from the next branching events by a number of generations of order n. Intuitively, 
the range of the competition kernel being smaller than the one of the mutations, the population 
fragments into small patches that have few interactions with each other. When a^ < au, the 
distinction of subfamilies in the genealogical tree is not so clear. The competition kernel has a 
large range and prevents the new mutants from escaping the competition created by the other 
individuals in the population. 

2.4.2 A variant of Adler's fattened goats: a spatial model 

In many models, trait or space play a similar role. Spatial models have been extensively studied 
as toy models for evolution (see Bolker and Pacala [SI [9] or Dieckmann and Law [IB]). Here 
we consider a spatial model where the competition exerted by past ancestors is softened. This 
model is a variant of Adler's fattened goats {e.g. [il|36]), but with an interaction in the drift 



term and not in the growth rate term). It corresponds to the choice of D 
with a > 0, r(t,y) = 1, b(t,y) = b and 



0, i^dids) 



d{t,y,X) 



Keiy'i 



yit)) 



JR'' 



K 



Xs{dy',dc')e-'''-^-'Us, 



(2.18) 



where, K^ is a symmetric smooth kernel with maximum at 0, for instance the density function 
of a centered Gaussian distribution with variance £. The death rate corresponds to the choice 
of U{t,y,y') = K^^yt — y'^). From the definition of the processes (see (j2.ip ). if y' belongs to the 
support of Xs{dy') then almost surely (a.s.) y' is a path stopped at s and for all t > s, y'^ = y'^. 
The goat-like particles consume resources at their location. When they arrive in a region 
previously grazed by the population, their death rate increases. The parameter a describes the 
speed at which the environment replenishes itself. The kernel K^ is the density function of a 
centered Gaussian distribution with variance e. The parameter K scales the carrying capacity 
and controls the mass of Xg . 


























1 ^^Jpj Wif 



Figure 3: Simulation for Adler's fattened goats, a = 10, e = 0.8, h = 0.75, a = 1. p = 1, K = 50 and 
n = 50. The 100 initial particles are started at location 1.5. 



Due to the form of the interaction, the goats spread quickly in the whole space and separate 



into families with very old MRCAs, as observed in Fig. [3] . The families become quickly disjoint 
and geographically isolated. 

3 The nonlinear historical superprocess limit 

We now investigate a diffusive limit for the sequence of processes X" defined by (|2.12|) . Let us 
firstly introduce a class of test functions which will be used to define the limiting process. 

Definition 3.1. For real C^-functions g on M+ x M and G on M respectively, we define the 
continuous function Gg on the path space Djjd by 



Gg{y) = G(^j g{s,ys)ds). (3.1) 



Let us remark that the class generated by finite linear combinations of such functions is stable 
under addition and separates the points, as proven in Lemma I A. II in Appendix |X1 Notice that 
if y is a cadlag path stopped at t G [0, T] then 



Gg{y) = G\^j g{s,ys)ds+ g{s,yt)dsy 



'0 Jt 

Also, in the sequel, the following quantity will appear for t € [0, T] and y € 

rT , rT 



V^Gg{t,y) = G'(^f g{s,ys)ds^ j A^g{s,yt)ds 

a J. f J. o 

g(s,y,)ds)^(y d,Ms,yt)ds) . (3.2) 



This quantity generalizes the Laplacian operator to an infinite-dimensional setting. For instance. 



if G{x) = X and if g{s, y) = g{y) does not depend on time, we get V'^Gg{t, y) = {T — t)Ag{yt). 



Note that Dawson ([l3], p. 203) and Etheridge ([H], p. 24) introduce another class of test 

functions of the form 

m 

'^iy) = Yl9j{yt,), (3.3) 

for m € N*, <ti < ■ ■ ■ < tm and Vj G {1, . . . , m}, gj G C^(M°', M). This class is not convenient 
when dealing with cadlag processes since the function ip is not continuous for the Skorohod 
topology. However, these test functions will be used when y is a continuous path. 
If y is a continuous path stopped at t then ip{y) = Y[T=i9j{ytjAt)- For a path y G C(M+,M'^), a 
time t > 0, we define 

m~l k m 

A^it,y) = Y, h.,t,^i[it){ll9jiytM n 9j)iyt)), (3.4) 

fc=0 i=l j=k+l 

where to = 0. 

The following lemma links the test functions ()3.ip and ()3.3|) . It is proved in Appendix Rl 



Lemma 3.2. Let ip he a test function of the form ()3.3p . Then there exists test functions (ipq) 
of the form (j3.ip such that for every y G C(M-|-,M°') and t G [0,T], the sequences {^q{y))qeN* 
and {T>'^ipq{t,y))q^^* are bounded uniformly in q, t and y and converge respectively to (p{y) and 
A^{t,y). 
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3.1 Main convergence result 

Let us assume that the mitial conditions converge: 

3Xo € AipO^Rd), hm Xq = Xq for the weak convergence. (3.5) 

The main theorem of this section states the convergence of the sequence (X")„gN*: 

Theorem 3.3. Assume (|3.5p and (|2.14p . Then the sequence {X^)n^^* converges in law in 
1D)(M+, A^i7'(D]gd)) to the superprocess X G C{M.+ ,J^F{0^d)) characterized as follows, for test 
functions Gg of the form i3.1\) : 

Mf« = {Xt,Gg) - {Xo,Gg) - f f (pr{s,y)^ V^Gg{s,y) 

Jo Jd^^ ^ ^ 

+ j{s,y,X')Gg{y)^Xs{dy)ds (3.6) 

is a square integrable martingale with quadratic variation: 

rt 
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{M''^)t= j I 2r{s,y)Gl{y)Xs{dy)ds, (3.7) 



where T>'^Gg{t,y) has been defined in ()3.2p and "f{t,y,X*') defines the growth rate of individuals 
y at time t in the population X: 

j{t,y,X') = b{t,y)-d{t,y,X'). (3.8) 

n 

For the proof of Theorem 13. 3t we proceed in a compactness-uniqueness manner. Firstly, we 
estabhsh the tightness of the sequence (X")„gN. (Section l3.2p then use Prohorov's theorem and 
identify the hmiting values as unique solution of the infinite-dimensional martingale problem 
(j3.6p , (j3.7p . The main difficulties in the proof are due to the interaction between the individual 
genealogies, implying the nonlinearity in the limit given by the term d{t,y,X^). 

3.2 Tightness of (X")„gN* 

In this subsection, we shall prove that: 

Proposition 3.4. The sequence (/:(X"))„gN* is tight on'P{O{R+,MF0^Rd))). 

To do that, we use the following criterion characterizing the uniform tightness of measure- valued 
cadlag processes by the compactness of the support of the measures and the uniform tightness 
of their masses. Moreover, a tricky coupling will allow us here to deal with the nonlinearity. 

Lemma 3.5. The sequence of laws o/ (X")„gN* is tight in'P{D{M.^,Aip(I])^d))) if 
(i) VT > 0, Ve > 0, ^K C ©ud compact, 

sup P(3t G [0,r], Xf (Kt) > e) < e, 
where K^ is the complement set of 

Kt = {y\y'- I y e i^, t G [o,r]} c B^^. (3.9) 

(a) V Gg of the form h3.1\) . the family of laws (((X", Gg)))n&i* is uniformly tight in 7^(018^). 
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Sketch of proof. Recall that Oj^d embedded with the Skorohod topology is a Polish space. By 
Lemma 7.6 in Dawson and Perkins |15j, the set Kt is compact in Pjjd. Let (i^fc,T)A;GN* be the 
family of compact sets such that 

sup P (3t € [0,r], X^K^kT) > e/2') < e/2', 
Then, the set 

fcGN* 

is then relatively compact in }Ap{p^d) by the Prohorov theorem, since it corresponds to a tight 
family of measures. We can then rewrite Point (i) of Lemma 13.51 in VT > 0, Ve > 0, 3.^ C 
A^i?(D]gd) relatively compact, 

supP(3te [o,r], X[' i^)<e. 
neN* 

Moreover the class of functions Gg separates the point and is closed under addition. Thus Points 
(i) and (ii) of Lemma [33] allow us to apply the tightness result of Jakubowski (Theorem 4.6 [28]) 
and ensure that the sequence of the laws of {X'^')n&\* is uniformly tight in P(]D)(M+, 7V4ir(]D)iRd))). 

■ 

Proof of Proposition \3.4\ We divide the proof into several steps. 



Step 1 Firstly, we prove Point (ii) of Lemma 13.51 Let T > and Gg be of the form (j3.ip . For 
every t £ [0, T] and every A > 

n\{X^,Gg)\ >A)< l|g|looSup„gM*E((^r,l)) ^ (3_10) 

which tends to when A tends to infinity thanks to (|2.15p . This proves the uniform tightness 
of the family of marginal laws (X", Gg) for n € N* and t fixed. Then, the Aldous and Rebolledo 
criteria {e.g. [29]) allow to prove the tightness of the process {X'^,Gg). For e > and r] > 0, 
let us show that there exist no G N* and 5 > such that for all n > uq and all stopping times 
Sn<T^< {Sn + 6) AT, 

P(|A^f«-^^f«|>7?)<e and P(KM"'^«)t„ - (M"'««)sJ > t?) < e (3.11) 

where A^'^^s denotes the finite variation process in the r.h.s. of (j2.16p . 

Let us begin with some estimates. We fix t G [0,T], /i € W^ and a path y S Bj^d stopped at t. 
Using the Taylor-Lagrange formula, there exists Oy^t,h £ (0, 1) such that 

Gg{y\t\yt_+h)-Gg{y) = G(^j g{s,{y\t\yt, +h)s)ds) - g(^ j g{s,ys)ds) 

= G'[j g{s,y,)ds)A{y,t,h) + hj"(^j g{s,ys)ds + ey^t,h^{y,t,h))K{y,t,hf (3.12) 

where 

A(y, t,h)= f (^g{s, yt_ + h) - g{s, yt_ )) ds (3.13) 

converges to zero when h tends to zero. Another use of the Taylor-Lagrange formula for the 
integrand in A{y,t, h), shows the existence a family rjy^t,h,s £ (0, 1) such that 

A{y,t,h)= / [h-\7r,g{s,yt_) + -^h [Hess g{s,yt_ +r]y^t,h,sh)] hjds. (3.14) 
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Using (|3.12|) and (|3.14p . we integrate Gg{y\t\yt_ + h) — Gg{y) with respect to K'^{dh). Since 
k'^{h) is the density of the Gaussian distribution of mean and covariance u^Id/n, integrals of 
odd powers and cross-products of the components of h vanish. Thus: 

/ ^ (Gg{y\t\yt_+h) - Gg{y)^K''{dh) = G' U g{s,ys)ds^^ A^g{s,yt_ + 7]t,y,h,sh)ds 

+ 5r^"(y 9{s,ys)ds + ey,t^hMy,t,h)^Y,(^J d^^gis,yt_)ds) +^ 

where C is a constant depending on G, g, a^ and p but not on n. Therefore 

^2, 



hm n 

n— >'+oo 



/ iGg{y\u\yu-.h)-Gg{y)\K^{dh)-%-^V'Gg{u,y) =0. (3.15) 



Noting that for G and (7 in C^, 'D^Gg is bounded from the definition (|3.ip . we obtain the following 
upper bound: 



-pa 



2 



n\A''T^r-^Z'\)<^ [R'^{\\V^Gg\\^ + l) + {B + D)\\G\\^) supE( sup (Xf,!)) 

+ ||G||oof/i^d[0,r] sup E( sup (Xj", 1)2)1. (3.16) 

n&i* te[0,T] J 

For the quadratic variation process, 

E(|(M"'G«)t„ - (M"'^«)5j) 



loo^ 



(2f + —-) sup E( sup (X,", 1)) + '''^^ ' ^ sup E( sup (X,", l)^)] . (3.17) 

V n / „gN* te[0,T] " nGN* iG[0,T] ^ 



We thus obtain ()3.1ip by applying the Markov inequality and using the moment estimates of 
Proposition 12.61 



Step 2 Let us now check that Point (i) of Lemma l3.5l is satisfied. We follow here ideas introduced 
by Dawson and Perkins [15j who proved the tightness of a system of independent historical 
branching Brownian particles. Here, we have interacting particles, which makes the proof much 
harder. 

Let T G M+ and e > and ET be a compact set of Djgd. We denote by K* = {y*\y ^ K} C B^d 
the set of the paths of K stopped at time t. Recall that Kt defined in ()3.9p is the set of the 
paths of K and of their left-limited paths stopped at any time before time T. Let us define the 
stopping time 

S"^ = inf{t G M+ I X'^iKfr) > e}. (3.18) 

From this definition, 

P(3t E [0,r], X'^iK^) >e)= P(5," < T). (3.19) 

Our purpose is to prove that it is possible to choose K and no such that sup„>„^ P(5" <T) < e. 
We decompose ()3.19p by considering the more tractable X'^[{K'^Y) and write 



'>(5," < T) =P( 5," < T, X^iK' Y) > 2 ) + ^(S'^ < ^' X^i{K' y) < - 

<-E{X-^{{K^y)) + P(5," < T, X-^{{K^y) < I) (3.20) 
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by using the Markov inequality. We will show in Steps 3 to 5 that there exists rj € (0, 1) such 
that for n large enough, 



S^ < T, X^iK' r) <-)< P(5," < r)(l - r?). 



Together with ()3.20p . this entails that 



'{S^ <T)< 



2E{XJ^{{K^)^)) 



erj 



In Step 6, we will also prove that the compact set K can be chosen such that 



E(X?((i^^)^)) < 



e^rj 



This will conclude the proof. 



(3.21) 



(3.22) 



(3.23) 



Step 3 Let us prove (|3.2ip . Heuristically, the event {5" < T, X^{{K'^Y) < §} means that more 
than half the trajectories that exited K before 5" have died at time T. On the set {S"" < T}, 
y = ^ K^ implies y ^ K . Indeed, if a path stopped at time 5"" does not belong to K, then 
it is also true when it is stopped at T > 5". Thus on {S^ < T}: 



n/r„,5" 



-S"- 



xUiK'r)>xU{y'^ ^K'^}) 



(3.24) 



traits 











s^ 



time 



Figure 4: The compact K is the region between the two lines. The paths which are drawn correspond to 
the support of X^. X^{{K'^Y) counts the trajectories that do not belong to K between and T: here 
we have yi, y2, 2/4, 2/5 and ye. The quantity XJj^{{y " ^ K ^ }) counts the trajectories, at time T , that 
don't belong to K between and S" . Here, we have yi, y^, j/5 and j/g ; although 2/2 does not belong to 
K between time and T, it belongs to K between and S*". To obtain the trajectories accounting for 
X1j^{{y^^ ^ K }); "w^ count the descendants of the 3 points at time S"" corresponding to trajectories 
y^e ^ K^^ i^g g^j^ g^gp check that relation p.24p is satisfied. 



Hence: 



'S^ < T, X^iK'^Y) < I) <p(5r < T, X^({/? ^ K'"}) < I 
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(3.25) 



Our purpose is to obtain an upper bound of the form (1 — rj) with r] € (0, 1) for the probability 
under the expectation in the r.h.s. of ()3.25p . This term is the probabihty that the population 
issued from particles which at time S" satisfy y^^ ^ K^" , has a size smaller than e/2. In view 
of (|3.25p . we will work on the set {S" < T} until the end of the proof. 

Using (j2.15p and Markov's inequality, it is possible, for any r/ > 0, to choose A^ > large enough 
such that 

sup {X^,l)>N\Tsp] <7]. (3.26) 

S^<s<T ' 

We need to introduce some coupling with independent trajectories. Let us define the process 
{ZJ^{dy))t£R_^ as follows. The independent particles of Z" are started at time S^ with the 
trajectories of X^„ such that {y^^ ^ K'^^} and the initial condition is 

Z^^idy) = l^s^^KS^X^^idy). (3.27) 

Their birth and death rates are nr{t,y) and nr(t,y) + D + UN. By a coupling argument, we 
have 

P(X?({/?^i^5?})<e/2|J-sn 

<P((Z?, 1) < ^; sup {X^, l)<N\FsA+^{ sup (X,", 1) > iV | J-^j 

^ ^ S^<s<T ^' ^SV-<s<T 

<1-P( inf (Z,M)>^|J-5n)+r?. 

If we can exhibit 77 > such that for n large enough 



(3.28) 



e(%<t P( inf (Z^, 1) > ^ I Ts?)) > 27jF{S^ < T), (3.29) 

then from ([3:25|) and dS^H]), the r.h.s. of ^^25\i is strictly smaher than (1 - 2r? + r])F{S^ <T) = 
(1 - r/)P(S'^ < T), and (fOTT) will be proved. 

Step 4 Let us prove ([3:29]) . Notice that on the set {S^ < T}, 

(Zg„, 1) = Xg„({/? ^ K'"}) > Xg„({/? ^ Ks^.}) > Xg„({/" ^ Kt}) = X^s^iKir) > e. 

By coupling arguments (using deletions of particles in the initial condition Zgn), and since we 
are considering a lower bound with an infimum in ()3.29p . we can consider without restriction 
that {Zgn, 1) = {[ne] + l)/n, where [x] denotes the integer part of x. 

Let us establish a diffusion approximation of {Z'^n, , 1) when n is large. We know that for any 
t > 0, the process 

(Zg„+„ 1) = (Zg„, 1)-{D + UN) I {Z^s^^^, l)ds + M,"'^ (3.30) 

where M"^'^ is a square integrable martingale such that for all s < t, 

2RJ\z-Si^^, l)du < (M"'^), - (M"'^), = ^* (zg„+„, 2t{S: + n, .) + ^^^)du 

<{2R + ^^^) [\z-s?_,,,l)du. (3.31) 

Proposition 12.61 and adaptations of ()3.16p and ()3.17p allow us to establish that the laws of 
{{Z^n, , 1), (M"'^).) are uniformly tight in ^(©(M-i-jR^)). As a consequence, there exists a 
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subsequence, denoted again by {{Zg„_^ , 1), (M"'''^),)„gN. for simplicity, that converges in distri- 
bution to a limit, say {Z, A) where Z and A are necessarily continuous. Let us define on the 
canonical space 

Aft = Zt-e+ I {D + UN)Zsds, (3.32) 

Jo 

and prove that it is a martingale. Let 0<si<---<S/fc<s<t and let (pi,- ■ ■ (j)k be bounded 
measurable functions on M_|_. We define 

^{Z) = M^s,) ■ ■ ■ M2s,)[Zt -Zs+ f {D + UN)Zuduy (3.33) 

From dSSQD, E(*((Zg„^ , 1))) = 0. Similarly to Prop. E6](h), we can prove from the SDE ([3:3(1 
that E(sup4gro j'i(^5n_,_4, 1)'^) < +00 for any T > 0. Then the sequence (^((Z^„^ ,l)))„gN* is 
uniformly integrable and by the continuity of ^, lim„_j._|_oo ]E(\I'((Z^„_|_ , 1)) = E(^(.E)). Then we 
deduce that K(^{Z)) = 0. Hence, Af is a continuous square integrable martingale, and Theorem 
6.1 p. 341 in Jacod and Shiryaev p5] implies that its quadratic variation process is (A/") = A. 
Moreover, using the Skorokhod representation theorem (see e.g. [7] Th. 25.6 p. 333), there exist 
a random sequence {Z^'',A^'')neN' and a random couple {Z,A) defined on the same probability 
space, distributed respectively as {{Z^n_^_ , 1), (M"'-^).)„gN. and {Z,A), and such that 

lim (Z",1") = (Z,^) a.s.. (3.34) 



n—^+oo 



Then, from (|3.3ip . we have a.s. that for all < s < t, 

2R / Zudu<At-As<2R Zu du. (3.35) 

J s J s 

This implies (see e.g. Rudin [381 Chapter 8]) that A is a.s. an absolutely continuous function 
and that there exists a random J^^-measurable function p{u) such that Vu G M+, 2R < p{u) < 2R 
and 

At = p{u)Zudu a.s.. (3.36) 

Jo 

Therefore, there exists a standard Brownian motion [Bt)t£'M,+ such that almost surely: 

M^ = Zt-e+ [ {D + UN)Zsds = I Jp{u)ZudBu. (3.37) 

Jo Jo 

Now that the diffusive limit for {Z^n, , 1) has been obtained, let us return to (j3.29p : 

P( inf (Zr,l)>^|-F5n)%<T = IP( inf (Zg„+„,1)>^|^5?)%<T 

>P( inf {y\s\Z^)s+u>^-) , , ]15?<T. (3.38) 

Notice that for all y and s, 



. -- ,,,„,_ ,,^_ J) =1P( inf {y\s\Z),+u>^i:)=^yJ inf ^„ > J 

n^+oo V«G[0,T] ' ' ^ ^ 2/ \mG[0,T] ' ' ' ^ 2/ ''' \uG[0,T] 2 



lim P( inf^Xy|^!^")«+-> n) =1P( mL(y|s|^)s+«>7T) =IP|/,4 inf^,^„>-), (3.39) 



where the notation P^^^ reminds that the distribution of Z depends on p which may itself depend 
on (y, s). We need some uniformity of the convergence in (|3.39p . with respect to y and s. 
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For C > and {z,r) G D x [0, T], we denote by B{{z,r),() the open ball centered at {z,r) with 
radius (. There exists C > small enough such that for all {y,s) £ B{{z,r)X) and Z £ ID, 
^Sk{{y\s\2), {z\r\Z)) < e/4 where dgk is the Skorokhod distance on ED (see Proposition IB . 1 1 in 
appendix). As a consequence, for this choice of C and all n € N*, 

,JPo^,(*l^"'-" > I) < F{jn,JMZ'U^ > 1)^ (3.40) 

Let .^ > be a small positive number. Since the sequence of laws of {Z"^, 1) is uniformly tight, 
there exists a compact set K^ of Bjgd such that for n large enough, P((Z", 1) ^ K^) < ^. Since 
K^ X [0,T] is compact, there exists a finite sequence {zi,ri)i<i<M with M = M(^) G N* such 
that 

A/(0 

K^x[0,T]c \J B{izi,ri),C). 

i=l 

With an argument similar to (|3.39|) . there exists for each i € {1, . . . , M}, an integer rii such that 
for all n > Ui, 

inf (zi\ri\Z"')r-+u>—)>-^zr( inf Z^ > —) . (3.41) 

Hence, thanks to ()3.40p . we obtain that for all {y, s) € K^ x [0, T] and n > maxx<j<jvf(^) n^, 

inf (y|s|^")s+„ > -) > min pf inf (zi|ri|^")r.+« > ^ 

> min -P^r-Y inf Zu > —) =: 2ri(e). (3.42) 
Then, from ()3.38p and (|3.42p . the left hand side of (|3.29p has the lower bound: 

>Efpf inf (y\s\Z'')s+u > -) Is^^rhz^- Dgk.) 

> 27^(0 ^{s:<T, {Z",l)eK^). 

The term P(S'" < T, (^", 1) G K^^ in the right hand side converges to P(5'" < T) when ^ tends 
to zero, and there exists ^o > small enough such that this term is larger than P(S'" < T)/2 
for every ^ < ^o- (In case P(<S'^ < T) = 0, the proof is done and this is also true). Thus, for 
< C < Co, the left hand side in ([3:i3]) is lower bounded by r/(C)P(5^ < T). This proves ([3:29]) 
provided r] is positive, which we shall establish in Step 5. 

Step 5 Let us prove that ij defined in ()3.42p is positive. Since it is a minimum over a finite 
number of terms, let us consider one of the latter. For this, we consider (z, r) € ORd and our 
purpose is to prove that 

P,,,( inf Zu>^) >0. 

V«G[0,T] 4 / 

For M > 0, let us define the stopping time ?m = inf{t >0,Z> M} and let us introduce 

r,/2 = inf{t>0, ii<|} (3.44) 
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such that Fz^r{^^^s€[o,T]^s > s/2) = Fz,r{Te/2 > T). Our purpose is to prove that the latter 
quantity is positive. We wih change the law of the process to be able to compare this quantity 
to the one defined with a Brownian motion. Let A > 0. Prom Ito's formula: 



gA2:tA5M =e^^ 



+ / {X^R-X{D + UN))Zse^'ds+ / XJ p{s)Zse^'dBs. 

Jo Jo 



Taking the expectation and choosing N sufficiently large (A^ > {XR — D)/U), we obtain that 
E(exp(A-EiA<;jvf)) ^ exp(Ae). From ()3.37p and since 2R < p{u) < 2R, we can classically show 
that E(sup(grQ J.1 Zf) < +00, from which we deduce that limM^+oo?M = +00 a.s. Moreover, 

it follows by Fatou's lemma that for any t G [0,T], E(exp(A2i)) < exp(Ae) and by Jensen's 
inequality and Fubini's theorem, we get 

^Uk 2p(s) ^-"'M < / E(e 2R ''')ds<e ^a < +00. (3.45) 

"^ u 

Novikov's criterion is satisfied and Girsanov's Theorem applied to (j3.37p tells us that under the 
probability M defined by 



= exp( - / ^== dBu-- / ^^ du], 3.46) 

J't V Jo V/o(n) 2 7o p{u) 



dF.^r 

Z is a Brownian motion started at e. Then we have 



(r,/2 < T) + M(t,/2 >T) = 1 



e 



^^{Zr^r,^;) = |M(r,/2 < T) +E"(Zt]1.,/,>t) = e. 

If M(rj/2 < r) = 1, this yields thus a contradiction since we would obtain e/2 = e for the second 
equation. Thus, M(t£/2 <T)<1 and F^,r(''"£/2 < ^) < 1- This shows that r^ > 0. 

Step 6 It now remains to prove (|3.23|) . We follow Dawson and Perkins [15, Lemma 7.3]. For 
n E N*, we can exhibit, by a coupling argument, a process X^ constituted of independent parti- 
cles with birth rate nr{t, y) + b{t, y) and death rate nr{t, y), started at Xq and which dominates 
X". In particular, forjT > and for any compact set K cB^a, E(X^{{K^Y)) < E{X^{{K^Y)). 
The tree underlying X" can be obtained by pruning a Yule tree with traits in M , where a par- 
ticle of lineage y at time t gives two offspring at rate 2nr(t,y) + b{t,y). One has lineage y 
and the other has lineage {y\t\y + h) where h is distributed following K"'(dh). Using Harris- 
Ulam-Neveu's notation to label the particles (see e.g. Dawson [13]), we denote by V^'" for 
a G X = Um=o{^' l}™"*"^ the lineage of the particle with label a. Particles are exchangeable and 
the common distribution of the process y"'° is the one of a pure jump process of W^, where 
the jumps occur at rate 2nr{t, y) + 6(t, y) and where the jump sizes are distributed according 
to the probability measure ^(5o(c?s) + ^K"'{dh). We denote by P" its distribution starting from 
X € M . It is standard to prove that the family of laws (P", n G N*, x € C) is tight as soon as 
C is a compact set of M*^. 

At each node of the Yule tree, an independent pruning is made: the offspring are kept with prob- 
ability {nr{t,y) + b{t,y))/{2nr{t,y) + b{t,y)) and are erased with probability nr{t,y)/{2nr{t,y) + 
b{t,y)). Let us denote by Vt the set of individuals alive at time t and write a ;^ i to say that 
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the individual a is a descendant of the individual i: 

j=l a>-i J \ i=l a)^i 

i=l J ''%'i 

The bound e is an upper bound of the mean population size at T that descends from a single 
initial individual, when the growth rate is bounded above by B. For each e > there exists a 
compact set C of M and a compact set K of ID)]gd such that 

sup XJ(C") < e and sup sup P;;^((K^)'^) < e, 

ngN* ngN* yo€C 

which concludes the proof. ■ 

3.3 Identification of the fimiting values 

Let us denote by X € C(M+, A^i?(D)) a limit point of {X^)n£n*- Our purpose here is to 
characterize the limiting value via the martingale problem that appears in Theorem 13.31 
Notice that the limiting process X is necessarily almost surely continuous as 

sup sup \{X,-,^)-{X^ ,^)\<-. (3.47) 



For the proof of Theorem 13.31 we will need the following Proposition, which establishes 
the uniqueness of the solution of (|3.6|) - (|3.7|) . Since the limiting process takes its values in 
C{[0,T],Mf{C{R+,M.'^))), we will use the test functions in ([331) instead of (IXB . 

Proposition 3.6. (i) The solutions of the martingale problem (j3.6p - (|3.7p also solve the following 
martingale problem, where ip is a test function of the form ()3.3p and A has been defined in ()3.4p : 



/o 
is a square integrable martingale with quadratic variation 



Mt = {Xt,^) - (Xo,^) - I I [pr{s,y)'^ ^^{s,y) + ^{s,y,X')^{y))X,{dy)ds (3.48) 



{M^)t= / 2r{s,y)^\y)Xs{dy)ds. (3.49) 

(a) There exists a unique solution to the martingale problem (j3.48p - (|3.49p . 
(Hi) There exists a unique solution to the martingale problem ()3.6p - ()3.7p . 

In the course of the proof, we will need the following lemma, whose proof uses standard argu- 
ments with r{t,y) depending on all the trajectory (see (|2.4p ). 

Lemma 3.7. Let us consider the following stochastic differential equation on M driven by a 
standard Brownian motion B: 

Yt = Yo+ [ y/a^pr{s,Y^)dBs. (3.50) 

Jo 

There exists a unique solution to ()3.50p . 
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Proof of Proposition [X5l It is clear that (iii) follows from (i) and (ii). 

Let us begin with the proof of (i). We consider a function (/? of the form (|3.3p and we assume 
without restriction that the functions gj^s are positive. Firstly, let us show that (j3.48p defines a 
martingale. Proceeding as in the proof of Theorem 5.6 in [23], for k£f^*,0<si<...Sk<s<t 
and (pi, . . . 4>k bounded measurable functions on AdpC^Rd), we define for X £ B(M_(_, Aipi^Rd)) 
the function 

^{X) = MXs,)...MXs,){{Xt,v)-{Xs,v)- f f {pr{u,y)^A^{y) 



+ ^{u, y, X")) V9(y)) Xuidy) du] . (3.51) 

Let us prove that E(\E'(X)) = 0. We consider, for q £ N*, the test functions (pq{y) = Gg^{y) 
with G{x) = exp(x) and gg{s,x) = YlT=i ^^s{9j{x))k'^{'tj ~ ■s) where k''{x) is the density of the 
Gaussian distribution with mean and variance 1/q. Prom (J3.6p - (j3.7p and (|2.15p . the process 

Mf^ = {Xt,^,) - {Xo,ip,) - j I (pr{s,y)'^ V\,{s,y) 

Jo Jd^, ^ ^ 

+ ^{s,y,X')^g{y))x,{dy)ds (3.52) 



is a square integrable martingale, hence uniformly integrable. The latter property together with 
Lemma [32] implies that (j)i{Xs^) . . . (pki^sJiMf" - Mf'') converges to ^(X) in L^ and that 
E(^(X)) =0. 

Let us show that the bracket of M'^ is given by (j3.49p . By a similar argument as previously, we 
firstly check that the process 



{Xt,^f - {,Xo,^f - / / 2r{s,y)^\y)Xs{dy)ds - / 2{Xs,^) / (pr{s,y)^Aip{y) 

Jo Jo Jo Jo ^ -^ 

+ -f{s,y,X'My))Xsidy)ds (3.53) 
is a martingale. In another way, using Ito's formula and (j3.48p . 

{Xt,vf - {Xo,vf - {M^)t - I 2{Xs,v) I (pr{s,y)^'K^{y) 

Jo Jo ^ ^ 

+ j{s,y,X'My))Xs{dy)ds (3.54) 
is a martingale. Comparing (|3.53p and (|3.54p . we obtain ()3.49p . 

Let us now prove (ii). Let P be a solution of the martingale problem (j3.6p - (|3.7p and let X be 
the canonical process on C(M+, A4i?(Digd)). We first use Dawson-Girsanov's theorem (see [HI 
Section 5], [21t. Theorem 2.3]) to get rid of the non-linearities. This theorem can be applied since 



^( [ [ -f\s,y,X')Xsidy) ds) < +oo. 
^ -'0 Jd„, ' 



^-"md 

Indeed, 7 is bounded and E( suptg[o^T](^t) 1)^ ) < +00 by taking ()2.15p to the limit. Hence, 
there exists a probability measure Q on C(M+, A^ir(D]gd)) equivalent to P such that under Q, 
and for every test function (p of the form (|3.3p . the process 

Mr = {Xt,v) - (Xo,^) - f f ^''^'f''' A^is,y)Xsidy)ds (3.55) 
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is a martingale with quadratic variation (j3.7p . Thus, if there is uniqueness of the probabihty 
measure Q which solves the martingale problem (|3.55p - ()3.7p we will deduce the uniqueness of 
the solution of the martingale problem (J3.6p - (j3.7p . 

Let us now prove that the Laplace transform of X under Q is uniquely characterized using the 
solution Y of the stochastic differential equation ()3.50p . We associate with Y its path-process 
W e C(M+,C(M+,M'^)) defined by 

Wt = {Ytf,s)s&L+- (3.56) 

The path-process W is not homogeneous but it is however a strong Markov process with semi- 
group defined for ah s < t and aU ip G Cb(C(]R+,M'='),M) by 

SsMv) =^'^{^{Wt)\Ws = y'). (3.57) 

Moreover, the infinitesimal generator A oiW aX, time t is defined for all (p as in (|3.3p by 



A^it, y) = ^r(t, y)Aip{t, y). (3.58) 

Then it can be shown that the log-Laplace functional oiXt under the probability Q, L(s, t, y, p) = 
E'Q(exp(-(Xi, ip))\Xs = 6y.) satisfies 

L(s,t,y,v5) = e-^-*'^(2'), (3.59) 

where Vs^t^{y) solves: 

VsMv) =^{v{Wt) - I ^r{u,Wu){VuMWu)fdu \ Ws = y') 

=SsMy) - 1 ^^Ss,u [r{u, .) {VuM-)) ') iy)du. (3.60) 

Adapting Theorem 12.3.1.1 of [131 P.207], there exists a unique solution to ()3.60p . Indeed, let 
V^ and V'^ be two solutions. From (j3.6Up . we see that for i S {1,2}, 

sup \Vs,My)\ < sup \ip{y)\ = \\ip\\oo- (3.61) 

s,t,y y 

We have 

IK>(2/)-K>(2/)l=|^/*^.,«(^(^,-)((K>(0)'-(v;>(.))'))(y)'in 

Ss,u{riu,.)2Moo\VlM-)-VlM.)\)iy)du 

<pcr^ll'/'l|oo^ / Ss,u{\Vltip - Vltip\j{y)du. 

Uniqueness follows from the Dynkin's generalized Gronwall inequality (see e.g. [13t Lemma 
4.3.1]). 

In conclusion, the Laplace transform of Xt is uniquely characterized for every f > by 
Ejf^(exp(—(Xi, (/?))) = exp(— (Xq, Vo,t</')). Thus, the one-marginal distributions of the martin- 
gale problem ()3.55p - (|3.7p are uniquely determined and thus, there exists a unique solution to 
(I335H - (13771) . ■ 

It is now time to turn to the 
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Proof of Theorem \3.3[ Let X G C{Rj^,MF{^Rd)) be a limiting process of the sequence {X"')n£n* 

and let us denote again by (X"')„gN* the subsequence that converges in law to X. Since the 

limiting process is continuous, the convergence holds in D(M_|_, A^i?(D]gd)) embedded with the 

Skorohod topology, but also with the uniform topology for all T > {e.g. [6]). 

The aim is to identify the martingale problem solved by the limiting value X. We will see that 

it satisfies ()3.6p - ()3.7p which admits a unique solution by Proposition 13.61 This will conclude the 

proof. 

Firstly, we show that (|3.6|) defines a martingale. For A: G N*, let < si < . . . s^ < s < t and let 

^1, ... ^fc be bounded measurable functions on Mpi^Rd). Let G G C^(M, M), 5 E C°'^(M+ xM'^', M) 

and Gg be functions as in (|3.1|) . We define for X G B(R+, ^^^^(IDjgd)) the function 

<^{X) = MXs,)...MXs,){{Xt,Gg)-{Xs,Gg)- [ [ (pr{u,y)^V^Gg{u,y) 

+ 7(n, y, X-)Gg (y)) X^idy) du] . (3.62) 
Let us prove that E($(X)) = 0. From ^M . 

=e(^MK) ■ ■ ■ MK){K'''' - Ms'""')) 

=E(<I>(X")) +e{^MK) ■ ■ ■ MKJi^n + Bn)), (3.63) 

where 

An= I I riu,y)\n( [ Gg{y\u\yu + h)K^idh)-Gg{y)) (3.64) 



-P^V'Ggiu,y)]x:idy)du 
Bn= I I b{u,y) [ (Gg{y\u\y^ + h)-Gg{y))K^idh)X:idy)du. 

As X is continuous, ^ is a.s. continuous at X. Moreover, from |$(X)| < C(sup^<^(Xs, 1) + 
suPs<t(^s, 1)^) and Prop. 12.61 we deduce that (<I>(X"))„gi^* is a uniformly integrable sequence 
such that 

lim E($(X")) =E($(X)). (3.65) 

Estimates (|2.15p and ()3.15p imply that An and B^ defined in (|3.64p converge in L^ to zero. Hence 
(|3.63p and (|3.65p entail the desired result. The computation of the bracket of the martingale 
Ms is standard (see for instance [23]). 
The proof of Theorem 13.31 is now complete. ■ 

4 Distributions of the genealogies 

An important question is to gather information on the lineages of individuals alive in the pop- 
ulation. 

Firstly, remark in Proposition 13.61 we obtain a martingale problem introduced by Perkins |36| 
Th. 5.1 p. 64] and leading to the following representation result: under Q, (Xt)jgiR^ has the same 
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distribution as the historical superprocess {Kt)teR+ which is the unique solution of 

Yt{y) = Yo{y) + [ y'a^pr{s,Y^{y))dys (4.1) 

Ko = Xo, {Kt,ip)= [ ip{Y{yf)Ht{dy) (4.2) 

where {Ht{dy))t£M.^ is under Q a historical Brownian superprocess (see [I5j), and for (/9 in a 
sufficiently large class of test functions, of the form p.3p for instance. The stochastic integral 
/q ^ya'^pr{s, Y'^{y))dys is the H-historical integral constructed in [36, Th. 2.12]. 

4.1 Lineages drawn at random 

For t > 0, Xt is a random measure on Djjd and its restriction to P([0, t],]^'^) correctly renormal- 
ized gives the distribution of the lineage of an individual chosen at random at time t. Let us 
define Utidy) = Xt{dy)/{Xt, 1) such that for any measurable test function ip on Bj^di 

{f.u^) = §^. (4.3) 

For instance, choosing ip(y) = tA{y) for a measurable subset A C Bjgd, we obtain the proportion 
of paths belonging to A under the random probability measure fif Studying such random 
probability measure remains unfortunately a difficult task and we will also consider its intensity 
probability measure E/xj defined for any test function ip as (E/^t, (f) = E((//i, 99)). This approach 
has been used in cases where the branching property holds (for instance in [2]). The hope is to 
identify the infinitesimal generator of Markov processes that will give us the average distribution 
of a path chosen at random. 

Proposition 4.1. For t > 0, a test function ip as in (|3.3|) and fit defined in 

pr{s,.)a^ 

/O 

+E(y" {{fis,lis,;X'M-)) - {f^s,^){fis,l{s,.,X')))ds) 

+Ie(^ jS-^{{fis,v){f^s,r{s,.)) - {f,s,r{s,.M-)))ds) (4.4) 

Remark that the moment equation (j4.4p cannot be closed due to the nonlinearity. 



E((/zt, if)) =E((/xo, ^)) + e( r (/i,, ^!l^:l^A<^(5, .))ds) 



Proof. We consider (j3.6|) - (|3.7|) and Ito's formula: 



JO^d {Xs,l) 

-E(f [ j^^j(s,y,XnX,{dy)ds) 

Mil f [ T^^Ms,y)Xs(.dy)ds-2 f [ .^—^2r{s,yMy)Xs{dy)ds]) 

^•^LJo Jo^a {Xs,i^r Jo JDj^d \Xs,i) J/ 

=E((/xo, (/9)) + e( y (/i„ ^pris, .)a^Aip{s, .) + 7(5, .,XsM-))ds^ 



/O ' Wo 

This ends the proof. 
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In (j4.4p . we recognize two covariance terms under the probability measure fj,s- 

Coy f,^{r{s,.), If) = {fis,r{s, .)ip{.)) - {fis,f){fJ's,r{s, .)). (4.5) 

The covariance terms (j4.5|) can be viewed as the generators of jump terms. For example, 

Cov^^(7(s,.,Xs),v3) = (/i^,7(s,.,Xs)) / (/ ifiz) ' ' ' ^ ^, -ip{y))^j.s{dy). 

Heuristically, the distribution of the path of a particle chosen at random at time t is as follows. 
We consider diffusive particles with generator j?r(s, .)o"^A. At rate 2(/is5^(S) •))/(^si 1)) the 
particles are resampled in the distribution [ig biased by the function ri^s^ .) where individuals 
with larger allometric coefficients ri^s^ .) are more likely to be chosen. At rate (;Us,7(s,y, X^)), 
particles are resampled in the distribution //^ biased by the function 7(s, ., X*) which gives more 
weight to particles having higher growth rate. 

4.2 Case of constant allometric function and growth rate 

A particular case is when the allometric function r{s,y) and the growth rate 'y{s,y,X^) are 
constant and equal to R and 7 respectively. In particular, a null growth rate 7 = models the 
"neutral" case. In that case, we will show that the process of interest corresponds to a Fleming- 
Viot process with genealogies similar to the one introduced in |24]. Notice that this bears 
similarities with the connections established between superprocesses and (classical) Fleming- 
Viot processes (e.g. Perkins [35]). 
Let us introduce the notation 

{lJ,t<S)Ht,x{t,;-)) = / xit,y,z)fJ-t{dy)fitidz). (4.6) 

This quantity can help to describe the most recent commun ancestor (MRCA) of two individuals 
chosen at random in the population. 

Proposition 4.2. Ifr{s,y) = R and ^{s,y,X) = 7, then 

(i)CoVf,^{r{s,.),ip) =0 and Cow^^(7(s, .,X^), (/?) = 0. 

(ii) The lineage distributions under ]E(/U() are Brownian motions. 

(Hi) Let x{y, z) be a function of two variables such that y i-> xilJj z) and z i-> xku^ z) «^e of the 

form ()3.3p . For test functions 4>{fJ.) = (/i (81 /x, x)/(/J, 1)^; the generator of X becomes 

rFVi^^^ pR(T^ / X X ~,o) \ '^R f f , ,X{dy) , X X . 

where A^^'>x{y, z) = A{y ^ xiv, z)) + A{z ^ xiv, z))- 

In the expression of the generator L , we recognize a resampling operator, where the resampling 
occurs nonlinearly at the rate 2R/{X, 1). 



Proof of Proposition \4^ The proof of (i) is obvious. To prove (ii), let us remark that for any 
test function if of the form (|3.3p , (|4.4p yields 



t 

2' 



{E{nt),ip) = {Eifio),ip)+ / {K{fis),pR<j'A^{s,.))ds. 

JQ 
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Therefore, under E(/it), the hneages have the same finite-dimensional distributions as Brownian 

motions with diffusion coefiicient pRa^ . 

Let us now prove (iii). Let ip and tp be two functions of the form ()3.3p . Using Ito's formula, 

+ / [^^{fJ's'^fJ's,A{(pi>){s,.)) +2R -^— fis{dy))ds, 

Jo ^ ^ Jm. \^s,i^) ^ 

(4.7) 



where M'^''^ is a square integrable martingale. 

For a function xiy,z) = J2k=i ^kfkiy)ipkiz) with K & N and Xk € M, we can generalize ()4.7p 

by noting that 



jf 



{fj't(S>fit,x) = ^>^k{^J't,y^k){^^t,i'k), 

k=l 

K K 

A^^^x(y, ^) = XI >'k^^kiy)^kiz) + X ^k^k{y)^ilJk{z). 



fe=i fc=i 

Since every function % as in (iii) can be approximated by functions of the previous form for the 
bounded pointwise topology, the proposition is proved. ■ 

4.3 Feynman-Kac formula in the case without interaction 

Following the approaches of [21 [22] in the case without interaction, it is more convenient to 
renormalize (j4.3p by E((Xf , 1)) instead of {Xt,l). Let us define the random measure ft for a 
measurable function Lp on Bj^d by: 

Contrarily to /Uj, the measure vt is not a probability measure, but its intensity measure Ef^ is. 

Proposition 4.3. Assume that pa"^ r{s,y) = 1 for the sake of simplicity. In the case without in- 
teraction, that is if^{s, y, X'^) = 7(s, y), then we have the following Feynman-Kac representation 
for any test function ip as in (j3.3p .• 

where {Wt)teR+ ^-^ (^ standard historical Brownian process, with values in C(M-|_,R°'). 

Proof. Let us denote by Mj the probability measure defined by the right hand side of (j4.9p . 
Using Ito's formula, we check that both (Ki't)t£R^ and (Mt)tg]R^ solve the following evolution 
equation in {mt)t£M.+ '- for any test function tp as in (j3.3p . 

{mt, 9?) = (mo, ^) + 2 ^"^*' ^^(^^ ■))'^^ + / ((™''*' '/'^(^' •)) " ("^s> V^) ("is, 7(s, •))) <^«' 
which admits a unique solution. ■ 
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In the general case for r{s,y), p and a^, the standard historical Brownian motion has to be 
replaced with the path-process W introduced in ()3.56p and associated with the diffusion (|3.50p . 
The Feynman-Kac formula, in the case where there is no interaction but still possible dependence 
of the dynamics on the past, allows for instance easier simulation schemes for approximating 
quantities that can express in terms of Ez^j. The right hand side of 14. 91 can be approximated by 
Monte-Carlo methods without branching mechanism. 

4.4 Feynman-Kac formula in the logistic case 

Let us now establish a Feynman-Kac formula for the logistic case without past or trait depen- 
dence. 

Proposition 4.4. Assume that r(s, y) = R with pa^R = 1 and 7(5, y, X) = a — rj{X, 1), where 
a, r] > 0. For any test function ip as in ()3.3p . 



where (Wi)tg]R_^ is a standard historical Brownian process and N is the solution of the following 
equation: 

Nt = (Xo, 1) + / {2R + aNs - vNs)ds + / V2RNsdBt, (4.11) 

Jo Jo 

B being a standard Brownian motion independent ofW. 

Proof. Let 99 be a test function of the form ()3.3p and g G C^(IR+,M). From ()3.48p . we get that 
the process iVj = {Xt, 1) satisfies the following martingale problem: 

Ml =Nt -No- I {aNs - r)N^)ds (4.12) 

is a square integrable martingale with quadratic variation {M^)t = L 2RNs ds. Taking the 
expectation leads to: 

E{Nt) =E{No) + / E{aNs - ??iVf)ds. (4.13) 

Jo 

Using Ito's formula with ()3.48p . ()4.12p . ()4.13p and taking the expectation, 
E{N,) I V E(N„) ) J„ I \E{N,y "'2 *^'/ 

2 \E{Ns) ^ ' 



ds. 
If we define the measure vt on C(]R4-,]R ) x M_|_ by 



{vt,(f0g} =E\ 



E{Nt) 
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then {i't)t<=M.+ solves the fohowing equation 

(Pi, ip(S)g) ={vq, ip(^g) + I (Us, -'Kip ®g + ip® {g{n){an - rjn^)) 

+ (p® g'{n){an - rjn^) + -(f ® 2Rg" {n)n + (/? 2Rg'\ 

— {i>s,ip (S> g){i^s,l (S> {an — r]n )) ds. (4-14) 

Let Yt = exp ( f^ {aNs — 7]Ng)ds) . Using Ito's formula with the processes {Wt)teR+ and {Nt)teR+ , 
we get 



k[^^^^^^ [E{\~A,iW. 



E(y,: 



+ E[^ip{W,){aNs-vN])Ys 



^E(^ip{Ws)g{Ns)Ys)E(^aNs - vN^) 



WAn'dVAY. . - „, \ /^ in(W..)n"(NAY. - \l 

ds. 



^ ^,,=(H^^(,;, ,„«,_,,,..)), ,(l«|p,«^, 



The measure fj on C(M+,M ) x M+ defined by 

(utj^^g) =IE 



ip{Wt)g{Nt)Yt 
nYt) 



is hence another solution of (|4.14|) . The uniqueness of the solution of this equation provides the 
announced result. ■ 
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A Properties of the test functions (13.11 ) and (13.31 ) 



We begin with a lemma that will be useful to link the classes of test functions (j3.ip and 
This lemma also shows that the class of test functions ()3.ip separates points. 

Lemma A.l. For q € N*, recall that we denote by k'^{u) the density of the Gaussian distribution 
with mean and variance 1/q. Let g € Cfe(M'^,M). For G{x) = x and gq{s,x) = k'^(t — s)g{x), 
we have for all y E Bjgd and all t € [0, T] at which y is continuous that: 

ii™ Gg (y) =g{yt). 

Proof. First notice that all the Gg^ are bounded by Hs'lloo- Let e > 0. Since y is continuous 
at t, so is (7 o y and there exists a > sufficiently small so that for every s G {t — a,t + a), 
\g{ys) — g{yt)\ < £/2. We can then choose q sufficiently large such that 

k'^{t - s)ds < 



\t-s\>a 4||g| 
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Then: 



GgM-9{yt)\= / k'i{t-s){g{ys)-g{yt))ds 



<2\\g\\^£ k'>{t-s)ds + - ki{t-s)ds<e 

J\s-t\>a 2 J\s-t\<a 



We are now in position to give the: 

Proof of Lemma \3.SX We can assume without restriction that the functions gj in the definition 
of if (|3.5p are positive. Let us define for y € C(]R_) 



ipq{y) = exp M ^\og gj{ys)k'i{tj - s)dsy 



By Leninia lA.il the term in the integral is bounded uniformly in q and y and converges when 
q tends to infinity to YlT=i log S'j (?/tj ) • ^s a consequence, for every y G C(M+,M'^), the sequence 
{iPq{y))n<=N* is bounded and converges to 

m m 

exp (^ Yl i°g 9j {yt, )) = n 5j '^y^j ) = 'p(y) 

when q tends to infinity. Moreover, 

yZ^og gj{ys)k'^itj - s)ds) 

f ^ A,(log5j)(yi)/c«(tj - s)ds + J2{ Y. d^.{^oggj){yt)k'i{tj - s)ds 
* i=i j=i ''* i=i 



When g tends to infinity, we have by Lemma I A. 1 



hm V\q{t,y) =\{g,{ytM ^ 



g-^+co 



^x9j{yt 



E 



{dx,9j{yt)Y 



i=i 



j I tj>t 



5i(yt) ■^ 5](yt) 



+ E( E 

4=1 j I tjX 



dx,9j{yt) 

9j{yt) 



9j{yt)9k{yt) 



n( \( sr ^^9j{yt) , ^v^ v;^ dx,9j{yt)dx,9k{yt) 
9j^y^.)\ L. gjy^) +^L. L. 



This concludes the proof. 



Uj I t,>t9j{yt 



i=i 



B Technical result on concatenated paths 

Let H be the set of increasing bijections from [0,T] to [0,T], where T > 0. We recall that the 
Skorokhod distance is defined for y, z € B by: 



dsk(y, z) = inf max { ||y o A - z\\oo, sup 
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log 



/ A(t)-A(.) x 



(B.l) 



In the sequel, we consider y, z € D and s,r € [0,T]. Without loss of generality, we can 
assume that s < r. 



Proposition B.l. If dgk{y,z) < e and if s and r are sufficiently close so that: 



T-s- 



T — r . 



< max < log - , log -^ \ < £■ 

Then for all if; € ID, dsk{{y\s\w), {z\r\w)) < 3e. 

In the proof, we will need the following change of time Xq a 7i: 

r / T — r \ 

Xo{u) = -ulu<s + [r + — (n - s) lu>s- 

s ~ \ I — s J 

The bijection Aq depends on r and s. For u and v € [0,T], we have: 

-Ao('u) - Ao(f) 



(B.2) 



log 



u — V 



< max <j log ( - j , log ^ rp _^ 



)}• 



(B.3) 



(B.4) 



The right hand side converges to when r/s converges to 1, and is upper bounded by e under 
the Assumptions ()B.2p of Proposition IB.li 

Lemma B.2. For all w G B. // ()B.2p is satisfied, then dsk{w o Aq, w) < e. 

Proof. The infimum in (jB.ip can be upper bounded by choosing A = Ag" , which is the inverse 

bijection of Aq: 

s / T — s \ 

Xq^{u) = -ulu<r + is+ rj. _ (u - r)] lu>r- (B.5) 

For such choice, we have: 



;{0,] 



dsk{w o Ao, w) < max <j 0, max ( log -, log — ])■ < e. 



;)} 



Let us now prove Proposition lB.lt 
Proof. By the triangular inequality: 

^Sk{{y\s\w),iz\r\w)) < A + B, where 

A = d{{y\s\w), {y o Ao|r|u; o Aq)) 
B = d{{y o Ao|r|u; o Aq), (2;|r|u;)). 

By Lemma IB. 21 A < e. For the second term, using Lemma lB.21 again: 

B <dsk(y o Ao, z) + d{w o Ao, w) < dsk(y o Ao, z) + e. 
Now, since dgk(y, z) < e, there exists X £ Ti such that ||y o A — zj|oo < 2e and 

X{u) - X{v) 



sup 

u,v<T 



log- 



u — V 



< 2e. 



Then, considering the change of time Aq o A: 

dsk(yo Ao,^;) < max i ||yo Ao o Ag"^ o A-2:||oo, sup log 

*■ u,v<T 

Xq o A(u) — Xq o X{v) 



(B.6) 



(B.7) 



Aq o A(n) — Aq o X{v) I 



< max< lly o A — zjloo, sup 

"- u,v<T 



< max ( e, 3e 



log- 



A(n) - X{v) 



+ sup 

u,v<T 



U — V 

Xiu) - X{v) 



log 



u — V 
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